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摘 Em. 了 解 植被 生长 对 气候 变化 的 响应 是 厘清 生态 系统 动态 关系 的 重点 。 基 于 1990 一 2018 年 气 
象 数据 和 归 一 化 植被 指数 (NDVI) ,应 用 偏 相关 分 析 与 地 理 探测 器 等 方法 ,分 析 了 在 生长 季 , 毛 乌 素 
沙 地 东南 缘 不 同类 型 植被 年 均 NDVI 的 变化 趋势 ,探讨 了 年 均 气 温 与 年 总 降水 量 对 各 类 型 植被 的 
影响 。 结 果 表 明 :(1) 1990—2018 年 生长 季 研 究 区 植被 年 均 NDVI 显 著 与 极 显著 增加 面积 达 


97.9% ,整体 生态 环境 质量 大 幅度 改善 。 


2005 年 之 前 植被 年 均 NDVI 增 速 缓 慢 , 此 后 以 0.011.a-: 的 


速率 发 生 了 突变 增加 ,其 中 灌 从 类 植被 年 均 NDVI 增 长 幅度 最 大 。(2) 2000 年 为 年 总 降水 量 与 年 均 
气温 的 趋势 突变 点 ,突变 前 年 总 降水 量 以 -5$.510 mm'.a-1 的 速率 减少 ,此 后 以 $.541 mm-a B 3 EE 
加 , 且 主 要 依赖 于 大 雨 雨量 的 增加 ;年 均 高 温 与 年 均 低温 在 突变 前 上 升 速率 分 别 为 0.122 Ca 与 
0.230 C-al, 突变 后 ,年 均 高 温 下 降 速 率 为 -0.014 Ca', 而 年 均 低 温 上 升 速率 为 0.022 7C a7, (3) 
在 植被 年 均 NDVI 缓 慢 增 长 阶段 (1990 一 2005 年 ), 年 均 低温 对 植被 影响 较 大 ,与 不 同类 型 植被 年 均 
NDVI 多 呈 最 著 正 相关 ;在 植被 年 均 NDVI 快 速 增 长 阶段 (2006 一 2018 年 ) ,年 总 降水 量 与 不 同类 型 
植被 年 均 NDVI 呈 显著 正 相 关 ,大 降雨 事件 的 频 发 使 得 降水 量 对 于 植被 的 生长 起 主导 作用 。 年 总 
降水 量 与 年 均 气 温 尤 其 是 年 均 低 温 的 交互 作用 是 促进 植被 生长 的 关键 。 

Xx 键 词 : 毛乌素 沙 地 ; 气候 因子 ; 归 一 化 植被 指数 (NDVI) ; 地 理 探测 器 ; 极点 对 称 模 态 分 解法 (ESMD ) 


文章 编号 : 


在 全 球 气候 变 暧 的 背景 下 ,各 区 域 的 水 循环 、 
于 湿 状况 发 生 了 显著 变化 ,这 对 世界 各 地 植被 的 生 
长 发 育 产 生 了 巨大 影响 "。 归 一 化 植被 指数 (NDVT) 
能 够 很 好 地 反映 出 植被 及 生态 系统 参数 的 变化 ”， 
且 对 气象 因素 的 变化 有 着 灵敏 的 指示 作用 ,被 广 
泛 应 用 于 植被 对 气候 因子 的 响应 研究 中 。 在 气候 
因子 中 ,降水 与 气温 对 植被 的 影响 最 为 突出 , Gao 
等 对 我 国 东部 区 域 研究 表明 ,植被 生物 量 主要 受 
降水 量 和 气温 的 影响 ,但 气温 对 植被 生长 的 作用 更 
强 ; 李 晴晴 等 ”发现 黄河 流域 植被 覆盖 度 与 气温 具 
有 显著 的 正 相 关 关 系 , 春 、 夏 、 冬 三 季 的 植被 覆盖 度 
与 气温 时 显著 正 相 关 , 与 降水 量 关 系 不 显著 ; 冯 颖 * 
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阐述 了 毛乌素 沙 地 植被 盖 度 与 降水 量 和 气温 均 呈 
正 相关 , 且 与 降水 量 相关 性 较 高 。 由 此 可 知 ,降水 
量 与 气温 对 植被 的 影响 作用 复杂 ,这 可 能 与 研究 区 
域 或 时 间 尺 度 有 关 , 也 可 能 是 未 对 人 研究 区 内 各 因素 
分 离 ,从 而 造成 了 结果 的 偏差 。 近 年 来 ,已 有 学 者 
从 不 同 角 度 说 明了 降水 与 气温 影响 植被 生长 的 机 
理 。 其 中 , 马 资 花 等 "提出 次 降雨 量 及 降雨 次 数 的 
分 布 共 同 决定 了 科尔沁 沙 地 一 年 生 植被 的 生物 量 
和 生物 多 样 性 ;高 漳 等 "发 现在 陕西 省 南部 ,日 最 
低 气温 的 升 高 能 够 促进 植被 生长 发 育 。 

毛乌素 沙 地 是 我 国 北方 重要 的 生态 屏障 ,其 水 
REZ ,生态 环境 脆弱 ,植被 动态 对 气候 变化 十 
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分 敏感 。 已 有 学 者 就 干旱 地 区 植被 NDVI 对 降水 量 
与 气温 的 响应 做 了 相关 研究 "2” ,但 其 多 从 降水 
量 均值 与 气温 均值 的 角度 加 以 分 析 ,未 能 进一步 分 
离 各 气候 因子 对 毛乌素 沙 地 植被 生长 的 作用 , 且 缺 
少 对 于 各 因素 影响 程度 的 量化 。 因 此 ,本文 基于 生 
长 季 植 被 遥感 数据 以 及 1:1000000 F Efi i d fe, 
利用 1990—2018 年 毛乌素 沙 地 东南 缘 内 5 个 气象 
站 点 的 的 气温 和 降水 数据 ,分 析 了 1990 一 2018 年 生 
长 季 内 年 均 高 温 、 年 均 低温 与 年 总 降水 量 的 变化 对 
不 同类 型 植被 NDVI 变 化 趋势 的 影响 , 旨 在 分 离 相 
关 因 素 ,阐明 气温 与 降水 量 的 变化 规律 及 其 对 植被 
生长 的 影响 机 理 , 并 利用 地 理 探测 器 模型 量化 各 因 
素 对 植被 生长 变化 的 作用 和 贡献 ,可 为 合理 制定 未 
来 生态 规划 与 反映 植被 对 于 气候 变化 的 响应 提供 
重要 依据 。 


1 研究 区 概况 


毛乌素 沙 地 东南 缘 (36°49’~39°27’N ,107°15'~ 
110955' E) , E H8 24 3.1x10* km ,行政 区 划 包 括 陕西 
省 榆林 市 定 边 县 、 靖 边 县 .横山 区 、 榆 阳 区 神木 市 5 
个 地 区 (图 1), 位 于 中 国 北方 农 牧 交错 带 。 该 区 域 
是 典型 的 生态 脆弱 区 ,更 是 干旱 ,大 风土 地 沙漠 化 
及 沙尘暴 频 发 的 多 灾 地 带 , 属 于 温带 半 干 旱 大 陆 性 
季风 和 气候 ,年 平均 气温 约 为 9.7 % ,年 平均 降水 量 约 
为 393 mm. 


2 数据 与 方法 


2.1 数据 来 源 

NDVI 数 据 来 自 GIMMS 3g v1 数据 集 以 及 中 
科学 院 资 源 环境 科学 数据 中 心 (http:/www.resdc.cn/ 
DOI) 发 布 的 中 国 年 度 植被 指数 (NDVI) 空 间 分 布 数 
据 集 。 其 中 ,CIMMS 3g v1 数据 集 时 间 跨 度 为 1981 
年 7 月 一 2015 年 12 月 ,空间 分 辨 率 为 8 km; 后 者 的 
时 间 跨 度 为 1998 一 2018 年 ,空间 分 辩 率 为 1 km。 对 
NDVI 数 据 的 预 处 理 包括 数据 格式 转换 ,通过 Arc- 
GIS 10.6 进 行 坐标 的 统一 .投影 以 及 裁剪 等 ,其 中 对 
于 分 辩 率 为 8 km 的 1990 一 1997 年 的 数据 集 ,通过 
重 采样 功能 与 掩 膜 裁剪 使 之 转化 为 1 km 的 数据 集 ， 
并 与 1998 一 2018 年 的 数据 集 相 整合 ,以 便于 后 期 
计算 。 

研究 中 使 用 的 气象 数据 为 中 国 气 象 数据 中 心 


人 


图 例 


。 气象 站 点 
0 80km 一 一 河流 
o — 县 区 分 界线 


植被 分 布 类 型 
mm 栽培 植被 mm 沼泽 上 DÉJA mem 草丛 


= 草原 
图 1 毛乌素 沙 地 东南 缘 植 被 分 布 


Fig. 1 Distribution of vegetation in the southeastern margin 
of Mu Us Sandy Land 
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的 地 面 气候 资料 ,1990 一 2018 年 的 日 值 数 据 集合 
(http://data.com.cn/site/index.html) , 以 分 析 年 均 高 
温 \ 年 均 低温 与 年 总 降水 量 在 各 时 间 尺 度 的 变化 趋 
势 及 其 对 NDVI 的 影响 作用 ,其 中 定义 年 均 高 ( 低 ) 
温 为 每 日 最 高 ( 低 ) 气 温 的 年 际 平均 值 ,定义 研究 区 
生长 季 为 每 年 5 一 9 月 。 按 照 国 家 气象 局 颁布 的 降 
水 强度 等 级 划分 标准 ,将 研究 区 24 h 降 水 总 量 分 为 
小 雨 (0.1~9.9 mm)、 中 雨 (10.0~24.9 mm) AA (25.0~ 
49.9 mm) 及 暴雨 (50.0 mm 及 以 上 ) ,为 了 评估 降雨 对 
植被 的 影响 ,本 次 只 保留 0.5 mm 以 上 的 有 效 降雨 "1。 

植被 类 型 数据 来 源 于 国家 自然 科学 基金 委员 
会 “中 国 西 部 环境 与 生态 科学 数据 中 心 ” 所 整理 的 
1:1000000 中 国 植被 图 ,毛乌素 沙 地 东南 缘 的 植被 类 
型 分 布 如 图 1 所 示 。 依 据 分 布 面积 从 大 到 小 依次 
为 :栽培 类 植被 (玉米 、 马 铃 莫 等 ) ,主要 分 布 于 定 边 
县 、 靖 边 县 ;草原 类 植被 ( 沙 党 、 从 生 和 不 草 等 ) ,主要 
分 布 于 榆 阳 区 、 神 木 市 ; 灌 丛 类 植被 ( 柠 条 Lb BU 
等 ) ,分 布 较 零散 ; 草 和 旬 类 植被 ( 茂 茂 草 等 ) ,主要 分 
布 于 定 边 县 北部 ; 阔 叶 林 类 植被 ( 杨 树 .柳树 、 榆 
树 ); 草 从 类 植被 ( 白 羊 草 ) ;这 漠 类 植被 ( 膜 果 麻 黄 、 
盐 节 木 等 ) 及 沼泽 类 植被 (芦苇 等 )。 

在 此 基础 上 ,通过 对 区 内 5 个 站 点 气象 数据 进 
行 克 里 金 插值 ,应 用 ArcGIS 10.6 提 取 每 种 植被 类 型 
所 在 区 域内 的 NDVI 气温 和 降水 数据 ,分 析 各 类 型 
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植被 与 气象 因素 的 关系 。 

2.2 研究 方法 

2.2.1 极点 对 称 模 态 分 解法 ”极点 对 称 模 态 分 解法 
(ESMD ) 在 气候 分 析 方面 具有 较 大 优势 "”。 本 文 利 
用 其 中 的 模 态 分 解 部 分 ,从 时 间 序 列 中 分 解 出 剩余 
变量 ,从 而 分 离 出 年 均 NDVI\ 年 均 高 温 、 年 均 低温 
与 年 总 降水 量 在 时 间 尺 度 上 变化 的 总 趋势 '“。 
2.2.2 Mann-Kendall 4 #4234 ”Mann-Kendall 趋 热 
检验 用 于 时 间 序 列 变化 趋势 分 析 '"”“。 设 时 间 序 列 
AN Xi Gz1, 2,…,n) ,F(x) 为 样本 XX 的 分 布 函 数 。 原 假 
WE Ho: F(x)2- FQ) ,为 n 个 独立 的 随机 变量 同 分 布 
的 样本 。 备 择 假设 Hi 为 双边 检验 。 对 于 所 有 i 中 的 
k jsn, Ake], XX MAY. BB at 
(5) 的 计算 公式 如 下 : 


n-l n 


s-YYXsew w 
E(S)=0 (2) 
Var(s)= Me Dem t (3) 
S-1. ($50) 
Var(S) 
Z=<0 (S=0) (4) 
—$*1 <9) 
Var(S) 
式 中 :5 为 检验 统计 量 ;Var(5) 为 方差 ;8(5) 为 均值 ; 
sgn 为 符号 函数 。 


在 双边 趋势 检验 中 ,对 于 给 定 的 置信 水 平 a ， 
若 |Z|> Z 。, 则 原 假设 Ho 不 可 接受 。 统 计 变量 Z 
2 


为 正 值 表示 增加 趋势 ,为 负 值 表示 减少 趋势 。|ZI 在 
大 于 等 于 1.96 或 2.58 时 ,表示 分 别 通过 了 信和 度 95% 
或 99% 显 著 性 检验 ,代表 结果 显著 或 极 显著 。 显 著 
性 程度 用 来 划分 研究 区 内 年 均 NDVI\ 年 均 气 温 与 
年 总 降水 量变 化 的 趋势 等 级 。 

2.2.3 趋势 分 析 法 ”趋势 分 析 法 用 于 反映 年 均 ND- 
VI\ 年 均 高 温 、 年 均 低温 与 年 总 降水 量 在 生长 季 内 
的 时 空 变化 趋势 中 ,计算 公式 如 下 : 


n n 
nx VixI,- Vixl, 
= i=l i=l 


04, e (5 ) 


式 中 :bw 为 各 像 元 年 均 NDVI、 年 均 高 温 、 年 均 低 温 
与 年 总 降水 量 的 变化 趋势 斜率 ;i 为 年 序号 ;n 为 总 
年 数 ;J; 为 第 i 年 的 年 均 NDVI .年 均 高 温 (%C) ,年 均 低 
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温 (%C ) 与 年 总 降水 量 (mm )。 

2.2.4 偏 相关 分 析 “ 偏 相 关 分 析 用 于 探索 植被 生长 
与 单一 气候 因子 之 间 的 联系 ,同时 剔除 其 他 气候 因 
子 的 影响 中 ,计算 公式 如 下 : 

站 Ls - T nx qua (6) 
(= rin Len 

式 中 :x 为 年 均 NDVI;y 为 年 均 高 温 (%C);z 为 年 均 低 
温 (%C);m 为 年 总 降水 量 (mm) ,rw 为 将 年 均 低 温 、 
年 总 降水 量 固定 后 年 均 高 温 与 年 均 NDVI 的 偏 相关 
系数 ,年 总 降水 量 年 均 低 温 与 年 均 NDVI 的 偏 相 关 
系数 同 理 。 若 之 0, 表 明年 均 高 温 .年 均 低 温 .年 总 
降水 量 与 年 均 NDVI 呈 正 相 关 ;r<0, 表 明年 均 高 温 、 
年 均 低 温 、 年 总 降水 量 与 年 均 NDVI 呈 负 相 关 ;r=0 
表明 年 均 高 温 .年 均 低 温 .年 总 降水 量 与 年 均 NDVI 
不 存在 线性 相关 性 。 

2.2.5 地 理 探 测 器 模型 ”地 理 探测 器 模型 用 于 探测 
空间 分 异性 ,是 揭示 自 变量 对 因 变 量 影响 机 制 的 一 
组 统计 学 方法 ”"。 本 文选 取 最 优 参 数 地 理 探测 器 
(OPGD)' 进 行 数据 处 理 , 在 Rstudio 软件 中 计算 每 
个 连续 型 因子 在 不 同 分 级 方式 .不同 间断 数 下 的 贡 
献 率 ,其 中 ,分 级 方式 根据 统计 规则 确定 (如 等 间 
PE, 自然 间距 .分 位 数 间距 .几何 间距 和 标准 偏差 间 
FE) ,间断 数量 设置 为 3~7 类 ,程序 通过 对 比 的 方式 ， 
选择 g 值 最 高 的 参数 组 合 (离散 方式 与 间断 数 )"*”， 
从 而 探测 影响 因子 对 NDVI 分 异 的 解释 力 。 

本 文 利用 地 理 探测 器 中 的 分 异 及 因子 探测 器 
评估 了 年 均 高 温 、 年 均 低温 与 年 总 降水 量 对 NDVI 
的 贡献 率 , 利 用 交互 作用 探测 器 分 析 了 年 总 降水 量 
与 年 均 高 温 \ 年 均 低 温 的 交互 作用 。 

分 异 及 因子 探测 器 :用 于 探测 不 同 植被 类 型 
NDVI 的 空间 分 异性 ,以 及 各 评价 指标 X 对 植被 ND- 
VI 变化 的 影响 力 大 小 ,用 g 值 度量 ,计算 公式 如 下 : 

Y No 

ei u 
式 中 ;有 =1, 2,…, 工 为 影响 因子 的 分 层 数 ;NN 与 N 分 
别 表示 在 1990 一 2018 年 生长 季 , 影 响 因子 年 均 高 温 
(Z)、 年 均 低 温 (三 ) 及 年 总 降水 量 ( 互 ) 在 产 层 及 整个 
研究 区 的 样本 数 ;o; 与 o 分 别 表 示 影 响 因 子 在 h 层 及 
整个 研究 区 域 NDVI 值 的 方差 。g 的 值 域 为 [0, 1], 
q70 表示 因子 与 NDVI 无 任何 关系 ;g=1 则 表示 因子 
可 以 完全 贡献 于 NDVI 的 空间 分 布 。 
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交互 作用 探测 器 :用 于 识别 评价 指标 之 间 的 交 
互 作用 , 即 评价 指标 承 与 无 共同 作用 是 否 会 增加 或 
减弱 对 植被 NDVI 的 解释 力 。 通 过 对 指标 各 自 g 值 、 
g(XNX2) 与 gi)+q(N%) 做 比较 ,将 交互 作用 类 型 分 为 5 
fh GR D", 


X1 交互 作用 类 型 
Tab.1 Interaction types 


交互 作用 类 型 4 值 关系 
非 线性 减弱 q(XNX,)<Min[q(X), q(X2)] 
单 因 子 非 线性 减弱 ”Min[g(X), g(X2)|<q(XiNX2)<Max[q(X,), qX) 
双 因 子 增强 qUXG X)» Max[q(X3), q(X:)] 
独立 XNXX) +X) 
EE ANRI X 5 XE HEH HI TERE 5 Min [ gX), qX) RA 
X, 5 X rp vA RE eME s Max [4 (X1), q(X2) RIR X X P BR oa 


度 最 大 值 ;g(Xi)+q(X;) 表 示 针 与 和 的 贡献 度 累 加 值 。X, 为 年 均 高 温 ; 
郊 为 年 均 低温 。 


3 结果 与 分 析 


3.1 气温 与 降水 的 变化 

时 间 上 ,在 1990 一 2018 年 生长 季 ,ESMD 结果 显 
示 2000 年 是 年 均 高 温 与 年 均 低 温 变化 的 转折 点 (图 
2a~b)。 其 中 ,在 1990 一 2000 年 ,年 均 高 温 的 增长 速 
率 为 0.122 Ca', 年 均 低温 的 增长 速率 为 0.230 C- 
a '; 在 2001 一 2018 年 ,年 均 高 温 的 下 降 速 率 为 
-0.014 Ca ,年 均 低温 的 增长 速率 为 0.022 C -a's 
对 于 不 同类 型 植被 所 在 区 域 而 言 ,ESMD 结果 显示 ， 
在 生长 季 , 各 类 型 植被 所 在 区 域 的 年 均 高 温 均 经 历 


在 区 域内 的 年 均 高 温 高 于 其 他 类 型 植被 (图 3a)。 
在 1990 一 2000 年 ,各 类 型 植被 所 在 区 域 的 年 均 低 温 
均 升 高 2 左右 ,2000 年 后 灌 从 与 栽培 类 植被 的 波 
动 幅度 小 于 0.8 C( 图 3b)。 

对 于 年 总 降水 量 而 言 ,ESMD 结果 显示 2000 年 
为 年 总 降水 量 趋势 的 突变 点 ,2000 年 之 前 ,年 总 降 
水 量 呈 下 降 趋 势 (-5.510 mm- a^) ,而 2000 年 后 年 总 
降水 量 显 著 增 加 (5.541 mm*a')( 图 2c)。 就 降雨 强 
度 而 言 ,在 2000 年 之 后 ,虽然 小 雨 雨量 所 占 平均 比 
重 (0.45) 最 大 ,但 其 占 比 呈 下 降 趋 热 (-0.003.a')， 
而 大 雨 雨量 占 比 却 呈 递增 趋势 (0.005.a) , 且 在 所 
有 降雨 强度 中 增长 幅度 最 大 ,这 表明 2000 年 后 年 总 
降水 量 的 增加 主要 依赖 于 大 降雨 事件 的 频 发 (图 
4)。 以 灌 从 类 植被 所 在 区 域 的 年 总 降水 量变 化 最 
明显 (图 3c)。 

空间 上 ,在 1990—2018 年 生长 季 ,研究 区 内 高 
温 多 年 均值 范围 为 26.33~27.62 ,变化 趋势 在 
-0.01~0.03 % .za 之 间 , 全 域 19.3% 的 面积 通过 了 显 
著 性 检验 ;低温 多 年 均值 范围 为 13.92~15.13 C , 增 
长 幅度 为 0.01~0.08 C- a`, 全域 79.7% 的 区 域 通过 
了 显著 性 检验 ;降水 多 年 均值 范围 为 386.12~511.17 
mm, 变 化 幅度 为 -0.72~4.95 mm: a” ,但 仅 有 3.96% 的 
面积 通过 了 显著 性 检验 ,说 明 研 究 区 的 降水 量变 化 
幅度 大 , 较 不 稳定 (图 5)。 
3.22 NDVI 的 变化 

各 类 型 植被 的 年 均 NDVI 均 呈现 显著 增长 趋 
势 , 其 中 灌 丛 类 植被 增长 幅度 较 大 (图 3d) ,为 
0.014*a (362), Hi ESMD 结果 可 知 ,研究 区 植被 的 


了 升 高 一 降低 一 升 高 的 变化 趋势 ,而 灌 从 类 植被 所 年 均 NDVI 在 2005 年 发 生 了 突变 增加 ,这 与 年 均 气 
(a) 年 均 高 温 (b) 年 均 低温 (c) 年 总 降水 量 
29 p —— o ESMD 趋 势 16 『 一 - 低温 均值 … ESMD 趋 势 600 [一 降水 量 均值 …- ESMD 趋 势 
一 一 阶段 趋势 一 二 阶段 趋势 一 一 阶段 趋势 一 二 阶段 趋势 一 一 阶段 趋势 一 二 阶段 趋势 
28 上 y- -0.014x426.793 500 F= 5.510x«313.26 
5 " 5 - = E 400 + 
上 pus 
i =0.022x+14.195 | A] 
m ái 2 m : TURO 2 Mal 
| y-0.122x26.039 il y=0.230x+12.132 200 + y-5.541x4298.94 
ag la R=0.333 — l m. R-0760 ， | R2-0.162 
1990 2000 2010 2020 1990 . 2000 2010 2020 1990 2000 2010 2020 
年 份 年 份 年 份 
注 :ESMD 为 极点 对 称 模 态 分 解法 。 下 同 。 
图 2 1990 一 2018 年 生长 季 年 均 高 温 \ 年 均 低温 与 年 总 降水 量 的 阶段 变化 趋势 
Fig. 2 Phase variation trends of annual extreme daily mean temperature, annual minimum daily mean temperature and total 


precipitation during the growing season from 1990 to 2018 
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一 一 栽培 植被 BR 一 一 草 旬 一 *- WA 
16 F (b) 年 均 低温 


28 F (a) 年 均 高 温 


27 
R 
mg 
26 
1 1 J 
1990 2000 2010 2020 
年 份 
500 F (o) 年 总 降水 量 
g 
Ë 400 
t 
x 
x 
300 
1990 2000 2010 2020 
年 份 


ge 
del 
me 
1990 2000 2010 2020 
年 份 
0.6 F (a) 年 均 NDVI 
z 
A 04+ 
z 
0.2 
1990 2000 2010 2020 
年 份 


注 :NDVI 为 归 一 化 植被 指数 。 下 同 。 
图 3 1990 一 2018 年 不 同类 型 植被 在 生长 季 的 年 均 NDVI 年 均 高 温 \ 年 均 低温 及 年 总 降水 量 的 ESMD 变化 趋势 


Fig. 3 ESMD variation trends of annual NDVI, annual extreme daily mean temperature, annual minimum daily mean temperature 


and total precipitation of different vegetation types during the growing season from 1990 to 2018 


—— 降水 强度 ”一 降水 天 数 一 - 降水 量 比重 


一 一 阶段 降水 量 比重 趋势 
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降水 量 比 重 


降水 量 比 重 


一 二 阶段 降水 量 比重 趋势 
(b) 中 十 
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图 4 1990 一 2018 年 生长 季 不 同 降水 强度 的 变化 趋势 


Fig.4 Variation trends of different precipitation intensities during the growing Season from 1990 to 2018 
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图 5 1990 一 2018 年 生长 季 年 均 高 温 \ 年 均 低温 与 年 总 降水 量 的 多 年 均值 .变化 趋势 与 显著 性 的 空间 分 布 


Fig. 5 Spatial distributions of multi year average, variation trend and significance of annual extreme daily mean temperature, 


annual minimum daily mean temperature and total precipitation during the growing season from 1990 to 2018 


温 及 年 总 降水 量 在 2000 年 的 突变 并 不 一 致 ,可 能 是 
因为 退耕 还 林 等 政策 的 实施 ,人 为 改善 了 研究 区 的 
植被 状况 ,造成 二 者 变化 的 不 同步 。 在 1990 一 
2005 年 生长 季 , 年 均 NDVI 增 长 幅度 为 0.001a , mi 
在 2005 以 后 ,NDVI 增 加 显著 ,增长 幅度 为 0.011.a7 
(图 6a)。 

在 空间 上 ,研究 区 内 NDVI 的 均值 分 布 范围 为 
0.14~0.64, 高 值 主要 分 布 在 定 边 县 东部 与 靖 边 县 西 
部 ; 低 值 主要 分 布 于 榆 阳 区 北部 .神木 市 西部 的 沙 
地 腹地 (图 6b)。 研 究 期 内 ,全域 年 均 NDVI 显 著 增 
加 与 极 显 著 增 加 的 区 域 面积 占 比 为 97.9% ,而 显著 
减少 与 极 显 著 减 少 的 区 域 面积 占 比 仅 为 0.1%。 


NDVI 减 少 区 域 主要 位 于 榆 阳 区 中 部 神木 市 西部 
与 靖 边 具 中 部 (图 6c)。 
3.3 NDVI 与 气温 和 降水 的 关系 

为 有 效 探 究 区 域内 植被 对 气温 与 降水 量 的 响 
应 规律 , 仅 对 占 比较 大 的 栽培 植被 .草原 、 灌 从 与 草 
馈 类 植被 进行 分 析 。 偏 相关 分 析 表 明 ,在 1990 一 
2018 年 生长 季 ,年 均 高 温 .年 均 低 温 与 年 总 降水 量 
对 各 类 植被 年 均 NDVI 的 影响 程度 存在 显著 差异 
( 表 3)。 控 制 年 均 低温 与 年 总 降水 量 的 条 件 下 ,年 
均 高 温 与 草原 . 灌 从 类 植被 的 年 均 NDVI 呈 显著 负 
相关 ;控制 年 均 高 温 与 年 总 降水 量 时 ,年 均 低温 与 
各 类 型 植被 的 年 均 NDVI 呈 极 显 著 正 相 关 ; 在 控制 
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表 2 1990 一 2018 年 不 同 植被 类 型 面积 占 比 与 其 NDVI 
在 生长 季 的 变化 趋势 
Tab. 2 Area proportion of different vegetation types and 
their change trends of NDVI during the growing season 
from 1990 to 2018 


年 均 高 温 与 年 均 低 温 的 条 件 下 ,年 总 降水 量 对 所 有 
类 型 植被 的 年 均 NDVI 均 显示 出 了 极 显著 的 正 相 
关 性 。 

不 同时 段 内 ,气象 因素 对 NDVI 的 影响 也 各 不 
相同 。 在 1990—2005 年 生长 季 , 年 均 低温 对 年 均 
NDVI 的 影响 最 显著 ( 表 3)。 其 中 ,在 控制 年 总 降水 
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量 与 年 均 高 温 的 条 件 下 ,年 均 低温 与 栽培 植被 . 草 
原 及 灌 从 的 年 均 NDVI 呈 显著 正 相 关 ; 在 控制 年 均 
低温 与 年 总 降水 量 的 条 件 下 ,年 均 高 温 与 灌 丛 的 年 
均 NDVI 呈 显著 负 相 关 ; 而 在 此 期 间 ,年 总 降水 量 与 
年 均 NDVI 并 没有 表现 出 显著 的 偏 相 关 性 。 但 在 


植被 类 型 ” 像 元 数 ” 占 比 /% ”变化 率 /a' 决定 系数 (R) PIE DE M ono 
栽培 植被 14465 48.3 0.010 0.864 «0.01 2008201 Je d 控制 ee 高 温 与 年 均 低温 的 
HEMA 678 23 0.014 0.840 <0.01 显著 正 相 关 。 同 时 ,控制 年 总 降水 量 与 年 均 高 温 的 
Hifa 407 13 0.009 0.96 ^ <0.01 条 件 下 ,虽然 年 均 低温 与 各 类 型 植被 的 年 均 NDVI 
LE $ M 0:006 pes x001 无 显著 正 相 关 , 但 对 栽培 植被 灌 从 和 草原 植被 的 
b 3 01 OON O84 — 400 生物 量 有 明显 的 促进 作用 。 这 表明 在 植被 NDVI 组 
荒漠 37 0.1 0.009 0.790 «0.01 慢 增长 阶段 (1990_2005 年 ), 年 均 低温 主导 了 植被 
沼泽 20 00 0.007 0.840 <0.01 Ren IS Bre mS APPT IRN i EF 


的 生长 ,而 在 植被 NDVI 快 速 增长 阶段 (2006 一 2018 
4E.) ,年 总 降水 量 对 年 均 NDVI 的 增加 贡献 最 大 ,并 
与 年 均 低温 共同 促进 各 类 植被 的 生长 。 

另外 ,地 理 探测 器 的 分 异 与 因子 探测 结果 显 
ZR ,在 1990 一 2018 年 生长 季 , 各 影响 因子 年 均 高 温 
(X) AES EE CX) CAE EK E (X) XT EX EM NDVI 
分 布 的 贡献 度 (g 值 ) 存 在 明显 差异 , 且 对 于 不 同类 


(a) NDVI 变 化 趋势 (b) NDVI EP py. (ONDVI 变 化 赵 势 显著 性 
03 p .. NDVISf& N ”显著 性 J 
---- ESMD 趋 势 图 例 
04 | 一 一 阶段 线性 趋势 
= ”| 一 二 阶段 线性 趋势 NDV 芍 值 
5 = 0.636 
Z o3 [770.0011:0.239 ^ 0.141 
R1-0.073 4 
ui à dd 765 -0.011x:0233 
0.2 + R1-0.847 
1990 2000 2010 2020 
年 份 


图 6 1990 一 2018 年 生长 季 时 间 和 区 域 尺 度 年 均 NDVI 变 化 趋势 


Fig. 6 Variation trends of annual mean NDVI on temporal and regional scales during the growing season from 1990 to 2018 


表 3 不 同类 型 植被 在 生长 季 年 均 NDVI 与 年 均 高 温 \ 年 均 低 温 及 年 总 降水 量 的 偏 相关 系数 
Tab.3 Partial correlation coefficients between annual NDVI and annual extreme daily mean temperature, annual minimum 


daily mean temperature and total precipitation during the growing season of different vegetation types 


年 均 高 温 年 均 低 温 年 总 降水 量 
植被 类 型 1990 一 。 1990 一 ”2006 一 1990 一 1990 一 2006 一 1990— 1990—  2006— 
2018 年 2005 年 2018 年 2018 年 2005 年 2018 年 2018 年 2005 年 2018 年 
栽培 植被 ” -0.33 -0.33 -0.36 0.66” 0.51' 0.52 0.53" 0.02 0.92" 
草原 -0.38' -0.42 -0.18 0.68" 0.58 0.35 0.50* -0.07 0.88" 
SH 草包 0.01 0.11 -0.17 0.55" 0.38 0.14 0.56" 0.31 0.81" 
HEA -0.49" -0.54 -0.31 0.67" 0.53: 0.44 0.54" -0.06 0.91" 


JÈ: NDVI 为 归 一 化 植被 指数 ;*、** 分 别 表示 在 0.05、0.01 水 平 ( 双 尾 ) 相 关 性 显著 。 下 同 。 
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型 植被 ,结果 也 各 不 相同 ( 表 4)。 整 体 上 ,年 总 降水 
量 的 贡献 率 最 大 ,其 次 为 年 均 低温 。 其 中 ,年 总 降 
水 量 对 草原 与 灌 从 的 生长 贡献 较 大 (g 值 分 别 为 
0.727 .0.622, P«0.05) ;年 均 低 温 对 草 多 类 植被 的 贡 
献 最 大 (g=0.713,P<0.01) ; 而 年 均 高 温 对 草 旬 类 植 
被 生长 贡献 最 大 ,但 不 显著 (q=0.390, P>0.05) 。 探 
测 器 交互 作用 结果 显示 ,年 均 高 温 .年 均 低温 与 年 
总 降水 量 的 交互 作用 对 植被 的 影响 均 为 双 因子 增 
强 , 结 果 均 大 于 单一 因素 对 于 植被 的 贡献 率 , ELTH 
较 于 年 均 高 温 而 言 ,年 均 低温 与 年 总 降水 量 的 交互 
作用 对 植被 的 影响 更 明显 ( 表 5)。 


表 4 单 因 子 影响 贡献 率 g 值 


Tab. 4 q value of single factor influence 


植被 类 型 影响 因子 4 值 Pi 
栽培 植被 年 均 高 温 /?C 0.287 0.341 
年 均 低温 /C 0.548 0.003* 
年 总 降水 量 /mm 0.541 0.023* 
草原 年 均 高 温 /人 0.354 0.396 
年 均 低温 /*C 0.535 0.062 
年 总 降水 量 /mm 0.727 0.000™ 
EX ESS eil C 0.390 0.310 
4FJ5 fil C 0.713 0.000* 
年 总 降水 量 /mm 0.399 0.599 
灌 从 年 均 高 温 /C 0.237 0.537 
年 均 低温 /C 0.501 0.007” 
年 总 降水 量 /mm 0.622 0.049' 


注 :4 为 因子 贡献 度 。 


表 5 评价 指标 交互 作用 结果 


Tab.5 Interaction results of evaluation indicators 


与 植被 的 生长 作用 关系 密切 。 其 中 ,比较 不 同 气象 
因素 与 各 类 型 植被 的 偏 相关 关系 及 其 贡献 度 可 以 
发 现 ,年 总 降水 量 对 毛乌素 沙 地 东南 缘 不 同类 型 植 
被 影响 更 大 ,并 且 降 水 量 越 丰富 越 能 够 促进 各 类 植 
被 的 生长 ,与 前 人 的 研究 结果 一 致 "2 。2000 年 
后 ,大 降雨 事件 的 频 发 导致 了 年 总 降水 量 的 变化 呈 
上 升 趋势 ,使 得 降水 量 对 于 植被 的 生长 起 主导 作 
用 ,这 证 实 了 Hesiler-Whiite 等 3\ 马 狗 花 等 "的 研究 
结果 , 即 在 干旱 地 区 ,增加 降雨 强度 有 助 于 植被 生 
产 力 的 提高 。 

在 1990 一 2018 年 ,年 均 低 温 与 不 同类 型 植被 年 
均 NDVI 呈 显著 正 相 关 , 且 在 1990 一 2000 年 年 总 降 
水 量 呈 下 降 趋势 的 情况 下 ,年 均 气 温 的 增长 ,尤其 是 
年 均 低 温 的 大 幅度 提升 , 仍 保证 了 植被 NDVI 的 增 
加 ,这 说 明了 低温 变化 对 于 植被 生长 的 重要 性 " 。 
在 温度 较 低 时 ,植物 存活 率 .叶片 伸 长 率 . 叶 绿 素 含 
量 和 地 上 净 初 级 生产 力也 最 低 习 ,同时 ,低温 胁迫 可 
大 大 限制 草本 植物 的 幼苗 出 苗 率 ,延迟 幼苗 物候 ”， 
而 且 冻 融 循 环 频率 的 增加 也 可 降低 多 年 生 植被 的 
存活 率 55 ,因此 ,在 中 高 纬度 区 的 毛乌素 沙 地 ,气温 
升 高 ,尤其 是 低温 的 升 高 可 改善 植被 的 生理 生态 表 
现 , 从 而 提升 植物 的 生物 量 。 

研究 区 内 灌 从 类 植被 年 均 NDVI 的 增长 速率 最 
快 ,其 所 在 区 域内 的 年 均 高 温 年 均 低温 及 年 总 降 
水 量变 幅 较 大 , 且 地 理 探测 器 结果 显示 年 总 降水 量 
与 年 均 高 温 或 年 均 低温 的 交互 作用 对 该 类 植被 的 影 
响 均 处 于 较 高 水 平 ,这 说 明治 从 类 植被 能 够 很 好 地 
适应 研究 区 内 气候 的 变化 并 保持 较 高 的 生物 量 ” ， 
建议 通过 人 工 修 复 措施 进一步 增加 沙 地 腹地 灌 从 


植被 类 型 交互 作用 交互 值 影响 力 x 3 、 
—— 一 一 一 - 类 植被 的 覆盖 度 ,从 而 使 毛乌素 沙 地 的 植被 恢复 持 
栽培 植被 Xfws 0.593 双 因 子 增 强 m 、 
Xf WX 0.802 双 因子 增强 续 向 好 发 展 。 
草原 Xx. 0.797 双 因子 增强 
A 
Xx. 0.921 双 因子 增强 5 结论 
254g] xx 0.587 双 因子 增强 
OH cs 本 文 基于 毛乌素 沙 地 东南 缘 5 个 气象 站 点 的 气 
灌 从 XX 0.634 双 因子 增强 温 、 降 水 与 NDVI 数 据 ,通过 ESMD 、 偏 相关 分 析 、 地 
xix. 0.827 双 因 子 增 强 理 探测 需 等 方法 前 明了 研究 区 内 1990—2018 年生 


注 :全 ,为 年 均 高 温 ;为 年 均 低 温 ;和 为 年 总 降水 量 。 


4 讨论 


本 文 利用 偏 相 关 性 分 析 与 地 理 探测 费 方 法 所 
得 结果 均 显 示 出 年 均 高 温 \ 年 均 低温 及 年 总 降水 量 


长 季 的 年 均 高 温 、 年 均 低 温 、 年 总 降水 量 与 年 均 
NDVI 的 变化 趋势 ,厘清 了 不 同类 型 植被 NDVI 对 各 
气象 因素 的 响应 规律 。 主 要 结论 如 下 : 

(1) 1990—2018 年 ,研究 区 植被 年 均 NDVI 显 
著 与 极 显 著 增 加 面积 达 97.9% ,而 年 均 NDVI 显 著 
与 极 显 著 减 小 面积 为 0.1% , 主要 位 于 榆 阳 区 中 部 、 
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神木 市 西部 与 靖 边 县 中 部 ,整体 生态 环境 质量 大 幅 
度 改善 。 在 2005 年 后 年 均 NDVI 以 0.011.a 的 速率 
发 生 了 突变 增加 ,其 中 灌 从 类 植被 年 均 NDVI 的 增 
长 幅度 最 大 。 

(2) 2000 年 为 年 总 降水 量 与 年 均 气温 趋势 的 突 
变 点 。 其 中 ,突变 前 年 总 降水 量 以 -5.510 mm-a WY 
速率 减少 ,此 后 以 5.541 mm*a'! 的 速率 增加 , 且 主 要 
依赖 于 大 雨 雨量 的 增加 ;年 均 高 温 与 年 均 低 温 在 突 
变 前 上 升 速率 分 别 为 0.122 C+ a^ 5j 0.230 Cal, R 
变 后 ,年 均 高 温 下 降 速率 为 -0.014 C- a, MEAR 
温 上 升 速率 为 0.022 Cea. 

(3) 分 离 气 候 因 子 对 植被 生长 的 相对 作用 ,能 
够 进一步 阐明 气候 变化 对 于 植被 生长 的 影响 。 各 
类 植被 生长 对 气候 变化 的 响应 较为 敏感 ,这 不 仅 说 
明了 水 热 条 件 会 影响 植被 活动 ,而 且 揭 示 了 其 对 植 
被 活动 的 影响 程度 。 其 中 ,年 均 低温 与 植被 的 生长 
多 呈 显 著 正 相关 ,在 植被 NDVI 绥 慢 增 长 阶段 
(1990 一 2005 年 ) 作 用 更 明显 。 年 总 降水 量 与 植被 
的 生长 呈 显 著 正 相关 ,在 植被 NDVI 快 速 增长 阶段 
(2006—2018 4E) ,影响 更 为 突出 。2000 年 后 ,大 降 
雨 事件 的 频 发 导致 了 年 降水 量 的 变化 呈 上 升 趋势 ， 
使 得 降水 对 于 植被 的 生长 起 主导 作用 。 在 整个 研 
究 期 内 ,年 总 降水 量 与 年 均 气温 ,尤其 是 年 均 低温 
的 交互 作用 是 促进 植被 生长 的 关键 。 
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Temporal and spatial variation of vegetation NDVI and its response to climatic 
factors in the southeastern margin of Mu Us Sandy Land 


HE Jungi, WEI Yan, GAO Wande, CHEN Yunfei, MA Yandong, LIU Xiuhua 


(School of Water and Environment, Chang’ an University, Key Laboratory of Subsurface Hydrology and Ecological Effect in Arid 
Region, Ministry of Education, Key Laboratory of Soil and Water Conservation and Ecological Restoration on the Loess Plateau of 


National Forestry and Grassland Administration, Xi’ an 710054, Shaanxi, China) 


Abstract: A central element in understanding the dynamic relationship between ecosystems is achieved by 
understanding the response of vegetation growth to climate change. Based on the meteorological data and 
normalized difference vegetation index (NDVI) from 1990 to 2018, the variation in the NDVI of different types 
of vegetation in the southeast margin of Mu Us Sandy Land of Shaanxi, China during the growing season was 
analyzed via use of a partial correlation analysis and geographic detector data, and the effects of temperature and 
precipitation on various types of vegetation are discussed. The results show that: (1) In the growing seasons from 
1990 to 2018, the annual NDVI for vegetation across the whole region increased significantly by 97.9%, and the 
overall eco-environmental quality was significantly improved. The growth rate of the annual NDVI of vegetation 
was slow prior to 2005, and then it increased sharply with a rate of increase of 0.011 .a '; across the observed 
increases, the annual NDVI related to shrub vegetation increased the most. (2) The year 2000 represented a 
transition year; a significant change in the behavior of the total precipitation and annual air temperature was 
observed in the year 2000. The total precipitation was seen to decrease at a rate of -5.510 mm - a`' before this 
transition year and increase with a rate of 5.541 mm-a ' after this transition year; these values were dominated by 
heavy rainfall events. Prior to this abrupt change in the year 2000, the annual rate of increase of annual extreme 
daily mean temperature and annual minimum daily mean temperature were 0.122 °C .an and 0.230 C » a', respectively. 
After the year 2000, the average annual decreasing rate of annual extreme daily mean temperature was -0.014 °C: 
a ', and the average annual rising rate of annual minimum daily mean temperature was 0.022 °C .a . (3) In the 
slow growth stage of the NDVI (1990—2005), the annual minimum daily mean temperature was positively 
correlated with the annual NDVI of different types of vegetation, and had the greatest influence; in the rapid 
growth stage of annual NDVI (2006—2018), the total precipitation was positively correlated with the annual 
NDVI of different types of vegetation, and the frequent occurrence of heavy rainfall events meant that the total 
precipitation had a dominant effect on the growth of vegetation. This work finds that the interaction between total 
precipitation and annual air temperature, especially the annual minimum daily mean temperature, is the key factor 
in the promotion of the growth of vegetation. 
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